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Abstract 

It is suggested that a large deep underocean neutrino detector, given the 
presence of significant numbers of neutrinos in the PeV range as predicted 
by various models of Active Galactic Nuclei, can make unique measurements 
of the properties of neutrinos. It will be possible to observe the existence 
of the v T , measure its mixing with other flavors, in fact test the mixing 
pattern for all three flavors based upon the mixing parameters suggested by 
the atmospheric and solar neutrino data, and measure the nu T cross section. 
The key signature is the charged current v T interaction, which produces a 
double cascade, one at either end of a minimum ionizing track. At a few 
PeV these cascades would be separated by roughly 100 m, and thus be easily 
resolvable in DUMAND and similar detectors. Future applications are precise 
neutrino astronomy and earth tomography. 

1 The Double Bang Signal 

In thinking about the consequences of r production in DUMAND it has 
become clear that we may find a spectacular signature for r events. This 
depends upon the existence of 10 15 eV neutrinos in adequate numbers, as 
are in fact predicted from AGNs0 for example. Since the r mass is about 
1.8 GeV, a r of 1.8 PeV and with cr of 91 //m|3j would fly roughly 90 m 
before self destructing. The interesting signals are from the charged current 
quark interactions of z/ T 's. The signature, as illustrated in Figure 0, is: 

1. a big hadronic shower from the initial v T interaction, 

2. a muon like r track, and then 

3. a second big particle cascade (usually 3 times larger). 




To give some scale to this, the ratios of detectable photons from these 
three segments are roughly 10 12 : 2 • 10 6 : 3 • 10 11 . 

The charged r will be hard to resolve from the bright light and not very 
different times from the cascades, but simply connecting the two cascades 
by the speed of light will suffice to make an unambiguous association. This 
appears to be a unique signature for real r production by z/ r 's, thus "discover- 
ing" the u T , and infering the mixing angles for neutrino oscillations. Finding 
even one of these events would have significant implications. 

The point is that in the general energy range of a few PeV there exists 
a powerful tool for searching for r mixing, over an unequalled parameter 
space, with unambiguous identification of the r. We know of no other way to 
make a v T appearance experiment with the cosmic rays, no way has been 
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proposed for an accelerator experiment except for the use of emulsions making 
observations of relatively large 5m 2 , and no way of detecting z/ T 's except 
statistically at proposed long baseline accelerator experiments. (Indeed the 
costs of these endeavors are similar if one includes the preparation of the 
beams at the accelerators.) 

In the following we explore the physics implications of the observations 
of the double bang events in a little more detail, discussing the kinematics, 
the sensitivity to two and three neutrino mixing, and potential backgrounds. 

2 Essentially Full Kinematics from Double Bang 
Events 

One gets to measure the total energy of the incident neutrino and nearly 
the full kinematics of the double bang events. The first cascade yields the 
energy transfered to the quark, and the second cascade gives the energy kept 
by the r; the sum give the total incoming neutrino energy, and the ratio 
of the first cascade energy to the total energy provides the y value. The 
cross sections and < y > are almost equal for v T and v T at this energy ||]. 
Observing the y distribution is a check on the observations, and departures 
from expectations could signal new physics. In calculating the v T flux, the 
measured y distribution will permit correction for the potentially unobserved 
events near y = (no initial cascade) and near y = 1 (initial cascade most 
of the energy and the tau decays too close for resolution). The near equality 
of the cross sections for particle and antiparticle permits the total flux to be 
calculated independently of the mix in the cosmic beam. 

The threshold energy for discriminating two cascades will be determined 
by having a r that flies far enough so that the two cascades can be distin- 
guished, and so that there are no "punch through" events. This distance will 
be of the order of some few times the cascade length (order 10 m) , and thus 
our threshold for r detection via this means would be, as suggested above, 
about a PeV. With the physics limitation of several tens of meters, there 
will also be a detector dependent limitation depending upon detector density 
and response. 

We note that the observation of the double bang events presents the 
opportunity to measure the PeV v T cross section via the angular distribution 
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in the lower hemispere caused by attenuation through the earth (~ 90% in 
this energy range). For future studies of earth tomography, the potential of 
this process is great, since it does not depend upon convolution over the y 
distribution and muon range, as is necessary to extract informatiom from the 
upcoming muon flux alone. 

Also, given the enormous light output of the cascades one would expect 
that timing from the detectors (at intermediate distances, since nearby de- 
tectors of present design will surely be saturated) would give excellent vertex 
resolution, and thus the initial neutrino direction to a precision of order on 
1°. In principle, of course, having both cascades and almost all energy "vis- 
ible" one can deduce the initial neutrino direction with arbitrary precision, 
perhaps making optical precision ultimately possible in some future neutrino 
telescope. 

3 Deducing the Neutrino Flux Flavor Con- 
tent 

From the ensemble of measurements with a DUMAND like array we will 
have: 

1. The r rate from double bang events gives the v T + u T flux. 

2. Measuring the UHE muon flux permits calculating the + flux. 

3. The W~ resonant event rate yields the v e flux at 6.4 PeV. 

4. Measuring the cascade rate (as a function of energy) gives the sum 
of neutral current interactions of all flavors of neutrinos and charged 
currents without visible //'s and r's, that is mostly i/ e 's. 

If we write r = (Jcc/&nc, the ratio of charged to neutral current cross 
sections, and we note that the cross sections are nearly flavor and charge 
independent at this energy, we can summarize the four observations abov as: 



Ni = N T + N f 



(1) 



N 2 = N ll + N t 



(2) 
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N 3 = N S 



(3) 



iV 4 = (N e + N s ) • (1 + r) + (N„ + N p ) + (N T + N f ) (4) 

Combining these four equations we can then extract 3 flavor fractions 
(f e = {N e + N E )/N totah f, = {N^Np)lN totaU ^df T = (N T + N f )/N total ), 
and the antiparticle to particle ratio, N s /N e . 

We expect nearly equal numbers of particles and anti-particles, with half 
as many z/ e 's as i^'s. In some scenarios there would be few z/ e 's. The v T flux 
from any source is expected to be very small in any case. Most predicted 
scenarios indicate adequate flight time for the 7r's and the //'s to decay at the 
source. In that circumstance there is only one free number in the initial flux 
ratios, the initial 7r + to 7r~ ratio. 

The quantities to be observed do have different energy dependent sen- 
sitivity, which will complicate matters. We would expect that systematic 
uncertainties (in factors such as the effective volume, energy calibration, 
etc.) will be important in interpreting results. Detector specific simulations 
are needed to make explicit evaluations of these measurement possibilities. 
Nonetheless, it seems to us to be practical to extract the total flavor content 
of the cosmic flux with reasonable precision from this suite of observations. 

4 Sensitivity to Neutrino Oscillations 

Our 5m 2 sensitivity (from L/E) is then fantastic, going down to the order of 
10 -i6 eV 2 ( the distance is out to the AGNs, ~ 100 MPc). To determine the 
two neutrino mixing angle sensitivity limit requires detector specific simula- 
tions. The limitation has to do with the AGN neutrino flux magnitude and 
effective volume for these events, and will probably be limited by statistics, 
at least in the near future. We guess it will be no better than 0.01 in the 
best of situations (as with most experiments). 

To discuss the (to be observed) fluxes in terms of neutrino oscillations we 
make a number of simplifying (but reasonable) assumptions. Explicitly, we 
assume that: 

1. the initial fluxes are in the proportion z/^:z/ e :z/ T ::2:l:0 (as generally ex- 
pected), 
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2. there are equal numbers of neutrinos and anti-neutrinos (although this 
is not crucial and can easily be dropped in actual analysis), 

3. all the 5m 2 are above 10~ 16 eV 2 , and that sin 2 (Sm 2 L/ '4E) all average 
to 1/2, 

4. matter effects at the production site are negligible (this is reasonable 
since N e - ~ N e + and baryon densities are low), and 

5. that there are no large matter effects in the path to earth. The latter 
is reasonable in the 5m range of interest (1CT 2 to 1CT 6 eV 2 ). 

With these assumptions, let us first consider the case where the low energy 
atmospheric neutrino anomaly is accounted for by neutrino oscillations. The 
oscillations could be either <-> v e or <-» v r with Sm 2 of 1CT 2 — 1CT 3 eV 2 
and sin 2 28 ranging from 0.5 to 1.0. 

In both cases the v^j v e ratio should be modified exactly as in the atmo- 
spheric case: given that the u^/u e is found to be 0.6 of the expected value 
of 2, we expect exactly the same thing from the distant cosmic sources of 
much higher energy, v^.v e = 1.2:1. In the <-> v e mixing case this translates 
into v^.v^.v^.-.X. 64:1. 36:0 (using the earlier normalization to 3 total). In the 
other case, of <-> v T mixing, we predict v^.v e :v T v\. 2:1:0. 8. It is easy to 
see that in the two flavor mixing case, v e j can never be greater than one. 
Even when the mixing is maximal on finds u^:iy e :iy T :: 1.5: 1.5:0 in the <-> v e 
mixing case, and i/^:i/ e :i/ T ::l:l:l in the <-> v T case. 

We can now ask how expectations will change if we include neutrino 
oscillation solutions to the solar neutrino deficit. We can consider two distinct 
possibilities: 

1. If the v e mixing with is with with either other species with a small 
angle, as in the small angle MSW regime, then the above results remain 
unaffected. 

2. If all three flavors mix substantially and the Sm 2, s are in the range of 
10~ 5 to 10~ 6 eV 2 (corresponding to the MSW "large angle" solution) 
or 10~ 10 eV 2 (corresponding to the "long wavelength" case). 

The parameter space for a combined fit to the atmospheric and solar 
neutrino data has been given for each of the above cases by Fogli, et aZ.[0], 
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and by Acker, et a/.[|J, respectively. We have evaluated the survival and 
transition probabilities for the whole range of allowed values of the three 
mixing angles, and the results are shown in Figure 0. This figure plots the 
fraction of muon neutrinos versus the fraction of electron neutrinos, and thus 
where each point specifies a fraction of tau neutrinos (it is analogous to the 
color triangle). The initial expected flux (u /M :i> e :i> T ::2:l:0) is at = 0.66 
and / e = 0.33. 

We further observe that: 

1. Amost all combinations of acceptable mixing angles result in saturation 
values to be observed with AGN neutrinos which lie between the lines 
fe + ffi = 0.88 and f e + / M = 0.66. Hence a substantial number of 
v T events are expected (0.34 > f T > 0.12) for all situations which solve 
the solar and atmospheric problems with neutrino oscillations. 

2. Since in the case of two flavor mixing it is impossible to obtain f e / / M > 
1, observation of the data falling above the diagonal f e = line is 
clear evidence for three flavor mixing. 

5 Backgrounds: Almost None 

We do not know of any particle with the ability to penetrate ~ 10 4 gm/cm 2 
of matter before decaying or interacting (100 m is about 100 strong inter- 
action lengths in water), except for leptons. Of the leptons, neutrinos have 
interaction lengths so great as to be negligibly likely to interact at close range 
(< 10~ 6 at this PeV energy). Electrons will immediately radiate. Muons 
are the most likely to cause confusion. A muon of 2 PeV energy has a mean 
energy loss rate of 500 GeV/m, and an effective radiation length of about 
1 m. Since on the order of 1/2 of the energy loss to the medium is continu- 
ous on the scale of a meter (it is a fuzzy number because of the continuous 
distribution of energy transfer fractions), the muon track light output will be 
roughly 1000 times minimum. Thus, even though muon radiation fluctuates 
greatly, we can anticipate that confusion with tau decays will be small. Once 
again, detector dependent simulations are required to numerically evaluate 
the confusion probability. Whatever that probability is, and we expect it 
to be low at 2 PeV, it decreases rapidly with energy (short range muon 
radiation goes up, and tau decay length increases). Finally, studies of the 
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Muon Fraction 

Figure 2: The fraction of muon neutrinos versus electron neutrinos, allowing 
for a fraction of tau neutrinos. The expected initial flux is at 0.66, 0.33. Full 
and equal mixing would result in 0.33, 0.33 (and 0.33 taus). The points repre- 
sent calculations of results for various solutions to the solar and atmospheric 
neutrino problems. 
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path length distribution can confirm the observation of tau decays, requiring 
consistency with the known tau lifetime. 

An other potential background might be due to downgoing muons. The 
downgoing muons are strongly peaked near the zenith, while the double bang 
events should be uniform in direction in the upper hemisphere, but heavily 
depleted from below the horizon by earth attenuation. A stronger constraint 
is that events with this energy are not expected from downgoing muons. 
Bremmstrahlung events of a few hundred GeV are plentiful, but the flux falls 
very fast with energy, and negligible numbers are expected above 100 TeV[f| 
(the radiation length for this energy muon is order of 20 m in water.). Thus 
even if one only had order of magnitude precision in cascade energy recon- 
struction (and we should have much better for near contained events such as 
we are discussing), one easily eliminates muons from the surface. 

6 Old Idea, New Combination and Impact 

The above is not entirely a new idea. Several authors have written about r 
signatures in the past|J. The new recognition is in coupling the uniqueness 
of the double bang signature, the encouraging flux predictions from AGNs, 
plus the modern hints about neutrino oscillations, and thus being able to 
make some claim as to the region we can probe in mixing space. Further 
we now recognize that observation of this novel class of interactions makes 
possible more precise measurements of the neutrino cross section and earth 
tomography than have been thought possible, because one can determine 
the neutrino total energy and direction with a precision limited only by 
the detector. For one thing, this potential observation does give motivation 
to filling in the volume, somewhat, of the hypothetical 1 km 3 array, and 
perhaps itself could justify construction of that experiment. Also, given that 
these events are near the acoustic detection threshold, one may contemplate 
hearing the double clicks from such events at km ranges and higher energies. 
(Bottom and surface reflected pulses might be thought to be a confusion 
factor, but in fact most often we would expect them to be useful as a tool 
for reconstructing range and orientation). 

We remind the reader that all of the above requires the existence of sub- 
stantial numbers of PeV neutrinos, which matter should be resolved in the 
next few years by the AMANDA, Baikal, DUMAND, and NESTOR exper- 
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iments now under construction. If those ultra high energy neutrinos are 
present in expected numbers, then we believe that the observation of the 
double bang events, along with other previously discussed interactions, will 
lead to important particle physics measurements which cannot be carried out 
in any other way on earth. 
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